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Abstract 

In models of coupled dark energy and dark matter the mass of the dark matter particle 
depends on the cosmological evolution of the dark energy field. In this note we exemplify 
in a simple model the effects of this mass variation on the relic abundance of cold dark 
matter. 



We still do not know the origin and composition of the cold dark matter (CDM) in the 
universe. Recent precision measurements of the cosmological microwave background and of the 
large scale structure of the universe put a strict bound on the abundance of non-baryonic CDM 

CP: 

Q CDM h 2 = 0.113 ±0.008, (1) 

resulting from a fit of several measurements combined in the framework of a ACDM model with 
a running spectral index. 

Most probably CDM is made of particles (even though there are alternatives where dark 
matter is the manifestation of a fluid with a non-standard equation of state, such as the Chap- 
lygin gas and quartessence models [2J) and the default candidate is the lightest supersymmetric 
particle (LSP) of supersymmetric (SUSY) extensions of the standard model of electroweak 
interactions, which is stable if R-parity is conserved. 

The calculation of the LSP abundance in the universe has reached a very sophisticated 
level. Computer codes are now publicly available that take into account all the several LSP 
annihilation and co-annihilation processes that enter in the evaluation of its abundance today 
The comparison of the results of the computer codes with recent observations are placing 
strong constraints in the parameters of SUSY (actually, minimal supergravity models) [3]. In 
fact, accuracies of the order of 10% among differents codes are being sought j^j. 

We also know that the universe is accelerating today. There should exist a form of dark 
energy, comprising roughly 70% of the energy density of the universe, responsible for its ac- 
celeration. The simplest possibility, still consistent with cosmological data, is a cosmological 
constant. One tantalizing problem that arises in these models is the so-called coincidence prob- 
lem: why dark energy starts to dominate the universe only at recent times? This has motivated 
the study of models in which dark energy is coupled to dark matter 1 . 

In these models of coupled dark energy, the mass of the dark matter particle depends on 
the dark energy field and therefore it varies on a cosmological time scale. In this letter we 
point out in a general way what consequences this effect may have in the calculations of the 
cosmological abundance of cold dark matter, illustrating them with a particular simple model. 

There are several different models of coupled dark energy, sometimes referred to as VAMPs 
(VAriable-Mass Particles) in the literature. The mass of the dark matter particles evolves 
according to some function of the dark energy field <ft, as, for example, a linear function of the 
field [3 IH1 13 EH] with a inverse power law dark energy potential or an exponential function 
[TT) IT2] IT3*1 ITU EI] with an exponential dark energy potential. For instance, if the dark 
matter is a fermion one could have an interaction like g{(j))mQXXi where the fermion mass 
m x = d( ( i ) ) m o i s a function of the dark energy field. Since the dark energy field is dynamically 
evolving with time, one would have m x = m x (0(a)) = m x (a). 

In order to have a rough idea of the possible magnitude of the mass- varying effect, we 
will follow a more phenomenological approach and assume that the dark matter particle mass 

1 There are also recent models that couples dark energy to neutrinos so that the energy density of neutrinos 
tracks the dark energy density Pj 
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evolves with the cosmological scale factor a as [TT] : 

m,(a) = mfet"^ (2) 

where is the particle mass today and a = In a. This results in the following equation for 
the evolution of the CDM energy density pcdm- 

Pcdm + 3Hp C DM - 5(a)HpcDM = (3) 

where H = a/a is the Hubble parameter. Conservation of the total stress-energy tensor them 
implies that the dark energy density should obey 

p DE + 3Hp DE (l + w DE ) + 5(a)HpcDM = (4) 

where Wde is the equation of state of the dark energy fluid. 

Many coupled dark energy models presents a so-called scaling solution where the dark matter 
and dark energy densities scale as simple functions of the scale factor. Following Majerotto et 
al. [T7j we impose: 

^ DE oc a^ (5) 
Pcdm 

This condition implies in the following expression for 5: 

5o 



n° DE + (i-n° DE )a-t 

where 



(6) 



So = -n DE {Z + aw DB ) (7) 

The free parameters of this model are taken to be 0, DE , wde and 5q. From a fit with the 
supernova la gold sample [T%] . Majerotto et al. obtain [T?] : 

n° DE = 0.62; w DE = -1.9] = -1.5. (8) 

which implies in £ = 8.1. 

With these values we can illustrate the effects of cold dark matter mass variation and we 
find that it saturates rather quickly for a < 0.1 at the value: 

= 1.34 (9) 

m x 

In order to estimate the consequences of the mass variation on the relic abundance of 
cold dark matter we have to solve the Boltzmann equation for the comoving number density 
Y = n/s, where n is the number density and s is the entropy density, as a function of the 
variable x = m(a)/T: 



x dY Ta 
Y E q dx H 



Y 
Ye~q 



(10) 
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where Yeq is the equilibrium distribution and Y a = Ueq^av) is the total annihilation width 
written in terms of the thermally averaged annihilation cross section times the relative velocity. 
The approximate solution of the Boltzmann equation gives the following relation for the CDM 
abundance today [TU] : 



n 



CDM 



oc K 



rn 



(0) 



(11) 

where Y^ is the solution of the Boltzmann equation for large x, when the number density 
is frozen after decoupling. For an s-wave annihilation process, which we will consider below, 
(ctav) = ctq and the solution can be well approximated by 



Yoo oc 



Xp 



(F) 



(12) 



( fp\ 

where xf is the freeze-out value of the variable x and m x is the CDM mass at freeze-out. We 
have kept only terms which have dependence on the CDM mass. In fact, in the following we 
will neglect a logarithmic mass dependence in xp to obtain 
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At this point, the only particle physics assumption was that of CDM s-wave annihilation. 
Therefore, naivelly one would think that there is a simple linear difference due to the possibility 
of mass-variation. 

However, the dependence can be potentially larger. As an example, we will consider a non- 
relativistic annihilation process XX ~ > // with a vector interaction with a coupling constant 
g mediated by a vector boson of mass My. In this simple case, away from the resonance (for 
My 3> Tn) we have 



Mxx -> //) 



9' 



mf) 2 (l + ^/3) 



(14) 



where /3f = y 1 — mj/ (m^) 2 . The important fact is the dependence on (m^) 2 . Hence one 
would get: 



(0) 



CDM 
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c. m. 

1 L CDM 



(15) 



where ^cdm i s ^ ne usua l result of the CDM abundance today with a constant mass particle. 

Therefore, in this simple illustrative model a difference of a factor of (1/1. 34) 3 = 0.41 in 
the cold relic energy density can be easily obtained. In these times of precision measurements 
and ruling out of parameter space for well defined models, the varying-mass phenomenom can 
be of importance. However, one should keep in mind that in realistic models there are several 
annihilation channels with different dependence on the CDM particle mass [20 . 

In summary, models that couple dark energy with dark matter predict mass-varying dark 
matter particles. Here we have exemplified the effect of mass variation in a simple model. Of 
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course any realistic model can be solved to find the function m(a) and verify the effect on the 
abundance. It would be interesting to implement the possibility of dark matter mass variation 
in realistic computational codes of relic abundances. 
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